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We study photon blockade and anti-bunching in the cavity of an optomechanical system in which the me¬ 
chanical resonator is coupled to a two-level system (TLS). In particular, we analyze the effects of the coupling 
strength (to the mechanical mode), transition frequency, and decay rate of TLS on the photon blockade. The 
statistical properties of the cavity field are affected by the TLS, because the TLS changes the energy-level struc¬ 
ture of the optomechanical system via dressed states formed by the TLS and the mechanical resonator. We find 
that the photon blockade and tunneling can be significantly changed by the transition frequency of the TLS and 
the coupling strength between the TLS and the mechanical resonator. Therefore, our study provides a method 
to tune the photon blockade and tunneling using a controllable TLS. 
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I. INTRODUCTION 

Cavity optomechanics has attracted extensive theoretical 
and experimental research activity in the last decade IMl- It 
ranges from testing fundamental aspects of quantum physics 
and gravity to applications in quantum enginee ring, qu antum 
measurements ll^ and weak-force detection lll5l - ll^ . For 
example, experiments ifisll have demonstrated the quantum 
ground state and single-phonon control in a mechanical res¬ 
onator, which is coupled to a superconducting TLS. It has also 
been shown that the mechanical resonator can be used for fre¬ 
quency conversion lfl9l - l^ . By controlling the frequency and 
the time intervals of a pumping field, nonclassical states of the 
mechanical motion can be prepared by carrier and sideband 
transition processes 

It is known that photon control can be realized in optome¬ 
chanical systems via an analogue of electromagnetic ally in¬ 
duced transparency (FIT), well-known in quantum optics. For 
instance, it has been found that FIT and photon scattering can 
be used to tune photon transmission in optomechanical sys¬ 
tems iIztI - I^ . a TLS coupled to the cavity field of an op¬ 
tomechanical system can affect the photon transmission and 
lead to nonclassical effects for the cavity field ||33 I - It^ . When 
the TLS is a controllable superconducting qubit, we find that 
the FIT window of the optomechanical system can be changed 
by the superconducting qubit, or, in other words, that the me¬ 
chanical resonator can affect the absorption and dispersion 
of the circuit QFD system j^ . Moreover, the mechani¬ 
cal resonator of an optomechanical system can also interact 
with a TLS which can affect the ground state cooling of 
the mechanical resonator S, the nonlinearity of the cavity 
field iHl], and so on. When the cavity field in such a hybrid 
system is driven by a strong classical field and a weak probe 
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field, the splitting of the phonon energy levels leads to two- 
color FIT windows which can be switched to a single 
one by adjusting the transition frequency of the TLS. 

Photon control in an optomechanical system can also be re¬ 
alized via photon blockade and tunneling, which result from 
the nonlinearity of the cavity field. Photon blockade prevents 
subsequent photons from resonantly entering the cavity, while 
the photon-induced tunneling increases the probability of sub¬ 
sequent photons entering the cavity. Thus, photon block¬ 
ade corresponds to a single-photon transition process, while 
photon tunneling corresponds to two-photon or multi-photon 
transition processes. If an optomechanical system coupled to 
a TLS via a mechanical resonator, both the mechanical res¬ 
onator and the TLS can induce a nonlinearity in the cavity 
field, and thus they can be used to realize photon blockade 
and tunneling. Photon blockade has been studied in various 
systems, e.g., cavity OL D ll43l-l^ . circuit QFD and 

optomechanical devices ll5^58]. In addition to the single¬ 
photon blockade, multiphoton blockade was also studied the¬ 
oretically (see, e.g., 115^^ and references therein) and even 
observed experimental!v 11^l63l - l^ . However, to our knowl¬ 
edge, there is no study on how to control photon blockade and 
tunneling. 

In this paper, we study a method to tune photon blockade 
and anti-bunching in an optomechanical system via a TLS 
which is coupled to the mechanical resonator of the optome¬ 
chanical device. In such a hybrid system El, the dressed 
states formed by the mechanical resonator and the TLS af¬ 
fect the photon and phonon blockade of the optomechani¬ 
cal system. It is known that the eigenstates of phonons in 
an optomechanical devices are described by displaced Lock 
states if^ due to the phonon-photon coupling via the radia¬ 
tion pressure. Therefore, the dressed states in the hybrid sys¬ 
tem should be more complicated, because they formed by the 
displaced phonon states of the mechanical resonator and the 
TLS IfiSt]. If the mechanical mode and the TLS are in 
the ultrastrong coupling regime, the rotating wave approxi- 
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FIG. 1: (Color online) (a) Schematic diagram of a hybrid structure 
consisting of a TLS coupled to the mechanical resonator of an op¬ 
tomechanical system. The TLS (within the black dashed circle) is 
denoted by a yellow dot inside the oscillating mirror represented by 
a black spring. Here, | J,) and | t) denote the ground state and ex¬ 
cited state, respectively, of the TLS. The parameters u>q, u>a, uJb, and 
LOd denote the frequencies of the TLS, cavity field, oscillating mir¬ 
ror, and driving field, respectively, (b) A schematic diagram for the 
couplings in the hybrid system with dissipation. The TLS is cou¬ 
pled to the mechanical resonator by the Rabi-type with the coupling 
strength g, the mechanical resonator is coupled to the cavity field 
with coupling strength x, and the cavity field is driven by an external 
field with amplitude fl. Here, a (a^) and b (b^) are the annihilation 
(creation) operators of the cavity mode and mechanical resonator, re¬ 
spectively, and Gx = (T+ + (T_. ja, Jb, and jq denote the decay 
rates of the cavity field, the mechanical resonator and the TLS, re¬ 
spectively. 


mation (RWA) doesn’t work, the Rabi type interaction should 
be considered. In particular, the effect of strong and ultra¬ 
strong c oupling on the photon blockade is analyzed in many 
systems 1691 - 1^ . 

The model studied here is a combination of the usual pro¬ 
totype optomechanical models. Hybrid systems composed 
of a TLS coupled to the cavity field of an optomechani¬ 
cal system have been studied widely (see, e.g., the recent 
Refs. ||33 I - It^ and references therein). Specifically, we 
consider a standard Hamiltonian for two interacting oscilla¬ 
tors (i.e., optical and mechanical resonators) in which the 
mechanical oscillator interacts also with a two-level system 
(TLS). It is worth noting that the model studied here is non¬ 
trivial because the couplings between its constituent subsys¬ 
tems are nonlinear. For example, the interaction between 
the two oscillators is proportional to the photon number and 
the position of a mechanical resonator. This nonlinear inter¬ 
actions can induce nonlinearity of the oscillators. For ex¬ 
ample, as will be shown below, the optical oscillator, due 
to its interaction with the mechanical oscillator, can be ef¬ 
fectively described by a Kerr-type nonlinearity. It is known 
that the standard Kerr nonlinearity can induce various non- 
classical effects Sill- These include self-squeezing iIt^ 
1^ . generation of two-component if^ Hoi and multi- 
component llsll [S^l Schrodinger cat states, and photon anti¬ 
bunching (if the nonlinearity is driven). The latter is a signa¬ 
ture of photon blockade (also referred to as optical state trun¬ 
cation) |l43l|43,ll6l], as also studied here. 

The creation of photons due to the mechanical resonator 
(i.e., oscillating mirror, which causes time-dependent varia¬ 
tions of the geometry of our mesoscopic optomechanical sys¬ 
tem) can be interpreted as a result of the dynamical Casimir 
effect (DCE), which is also known as non-stationary Casimir 
effect or motion-induced radiation (from a dynamically de¬ 
forming mirror). As explained in Ref. dHt]: “The term ‘dy¬ 
namical Casimir effect’ is used nowadays for a rather wide 
group of phenomena whose common feature is the creation 
of quanta (photons) from the initial vacuum (or some other) 
state of some field (electromagnetic field in the majority of 
cases) due to time-dependent variations of the geometry (di¬ 
mensions) or material properties (e.g., the dielectric constant 
or conductivity) of some macroscopic system.” Specifically, 
we can interpret the occurrence of photon blockade in the 
studied system as follows: As mentioned above, the nonlinear 
interaction between the mechanical and optical resonators of 
our system can induce an effective Kerr-type nonlinearity of 
the optical resonator. This driven Kerr nonlinearity can result 
in photon blockade. Note that this driving is applied directly 
via the coupling of the mechanical and optical resonators (be¬ 
ing related to the DCE) and indirectly via the coupling of the 
mechanical resonator with the TLS. 

The DCE was studied in analogous systems in a number of 
recent works (see, e.g., Refs. lIssUS^ ). In particular. Ref. ll^ 
analyzed strong modifications of the cavity field statistics in 
the DCE due to the interaction with TLSs. Here, we study 
the effect of a single TLS on photon blockade. The light gen¬ 
erated via the DCE can exhibit various nonclassical proper¬ 
ties ll^ including squeezing llssl - f^ 1^l9ll] . 
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The paper is organized as follows; In Sec. II, we describe 
the theoretical model. In Sec. Ill, we write down the mas¬ 
ter equation and derive the analytical solution in the weak- 
pumping limit. The photon blockade is analyzed via the 
second-order degree of coherence in Sec. IV. We finally sum¬ 
marize our results in Sec. V. 


II. ENERGY LEVEL STRUCTURE OE THE HYBRID 
SYSTEM 

A. Theoretical model 

As schematically shown in Fig. [T] we study a hybrid sys¬ 
tem which consists of an optomechanical cavity coupled to 
a TLS with its mechanical mode. We assume that there is no 
direct coupling between the TLS and the cavity field of the op¬ 
tomechanical part. In this case, the Hamiltonian of the hybrid 
system can be written as 




hcJaCL^CL + huJbb% + —LJqaz 
+hg +h)ax- 


hxa) a (b^ + b) 

( 1 ) 


Here, a (a^) and b (b^) are the annihilation (creation) operators 
of the cavity field and the mechanical resonator, respectively. 
The frequencies of the cavity field and the mechanical res¬ 
onator are denoted by oJa and uJb, respectively. The transition 
frequency of the TLS is ujq. The Pauli operators az and ax 
are used to describe the TLS with the ladder operators a± de¬ 
fined by ax = a-\. + a-. The coupling strength between the 
cavity field and the mechanical resonator is x, and the param¬ 
eter g describes the coupling strength between the mechanical 
resonator and the TLS. 

It has been shown (e.g., in Ref. that the mechanical 
resonator can mediate a Kerr nonlinear interaction between 
photons of the cavity field in an optomechanical system. Thus, 
when the mechanical resonator is coupled to the TLS, the me¬ 
chanical resonator will induce the effect of the TLS on the 
nonlinearity of the cavity field. To see this clearly, we apply 
a unitary transformation, U = exp [—xa^a{h^ — 6)/a;;,], to 
Eq. O. Then the total Hamiltonian in Eq. O becomes 



Photon number 


FIG. 2: (Color online) Schematic diagram of the energy levels of the 
hybrid system when the TLS resonantly interacts with the mechani¬ 
cal resonator (the ground state energy of the TLS is assumed to he 0). 
Here |n) (n = 0, 1, 2, • • •) represent the Fock states of the photons. 
If the photon number is zero, the dressed states of the phonon and 
the TLS are: |1±) = \l,i) ± |0,t) and |2±) = |2,;) ± |l,t). 
Here | f) and | f) denote the eigenstates of the TLS which is 
not coupled to mechanical resonator. When the photon number 
is nonzero, the expressions of the wave functions for the dressed 
states are given by |m±(n)) in Eq. (|5j (n, m = 0,1, 2, • • •). and 
the corresponding eigenvalues are given in Eq. (0 (in the stable 
regime). Here (5i = A'(l) — Aq and 82 = A'(2) — 4Ao, where 
A'(n) = 2gx{ax)/uJb and Ao = 'x /i^b- 

Aq = loJb as the photon nonlinearity induced by the me¬ 
chanical resonator in the optomechanical system, then the to¬ 
tal photon energy levels shifts for the one-photon and two- 
photon states of the hybrid system are (5i = A'(l) — Aq, and 
82 = A'(2) — 4Ao, respectively, as schematically shown in 
Fig.m Here, both A'ItA = 2gx{o'x) j^b and {a^ depend on 
the photon number n . 


B. Eigenvalues and eigenstates 


Hq = huaO'^O' + h-^^axtt^ a + —Wqaz + hujbb^b 
UJb 2 

—h—a'aa'a + Hg (b^ + b) ax- 
UJb 


( 2 ) 


Besides the energy level shift —v^-y^jujb (n = 0,1, 2, • • •) 
induced by the mechanical resonator with the photon num¬ 
ber n, the interaction between the TLS and the cavity field 
through h2gx<Jxa^a/uJb also leads to a photon energy level 
shift 2gx{o'x)/uJb, which is twice the result found with the 
RWA 114^ . This interaction induces a new nonlinearity of the 
cavity field. For example, in the case of large detuning be¬ 
tween the cavity field and the TLS, the TLS can induce an¬ 
other photon-photon Kerr interaction term . If we define 


We now analyze the eigenvalues and eigenstates of the sys¬ 
tem when the cavity field of the optomechanical system is in 
a Fock state |n) with the photon number n. In this case, the 
quantity xn can be considered as an effective driving field for 
the coupled system of the mechanical resonator and the TLS, 
and the effective Hamiltonian of the mechanical resonator and 
the TLS for the photon Fock state \n) can be given, from 
Eq. ([T]i, as 

Hb = huJbb'^b + ^ujqaz + hg {b^ + b) ax - hxn {b'^ + b) ,(3) 

where the constant term huan has been neglected. 

Let us first study the eigenvalues and eigenstates when the 
mechanical resonator and the TLS satisfy the resonant inter- 
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action condition in Eq. (IHi, i.e. = ojb — ujq = 0- Under the 
RWA, the dressed state energy levels in the interaction picture 
can be given as 


E, 


n,m,± 




(4) 


Here, m {m = 1,2,...) denotes the phonon number. If the 
effective driving field is not very strong, that is, 2xn < g, 
the dressed states will be stable. Otherwise the phonons have 
large chances to transit to hi gh e nergy levels and the dressed 
states will be unstable l6^ . The energy levels of the 
dressed states are also functions of the photon number. The 
dressed states are always stable for the 0 photon state, how¬ 
ever, the dressed states become unstable when n > 3/(2x)- 
We can see from Eq. (HI that each energy level has an extra 
term compared with ±hy/mg of the common dressed states 
in the resonant interaction between the TLS and the mechani¬ 
cal resonator. Here, the splitting width of the dressed states is 
affected by the quantum states of both photons and phonons. 

If the photon number is zero, the eigenvalues correspond¬ 
ing to the dressed states in Eq. © become zLtl^/mg, and 
the corresponding dressed states can be written as \m±) = 
[|to, ),) ± |to — 1, t)] /which are the common dressed 
states of the resonant interaction between the TLS and the 
mechanical resonator, as schematically shown in Eig. |2] The 
dressed states, corresponding to the eigenvalues in Eq. ©, can 
be written as ibl, 

|m±(n)) = l)\P) 

±i\r],/3{En,m,±);m)\M)]. (5) 


In circuit QED, the standard photon blockade is signifi¬ 
cantly changed by the ultrastrong coupling between the cav¬ 
ity field and the TLS lHH]. Since the sideband-transition pro¬ 
cesses in optomechanics usually accompany the absorption or 
emission of phonons, the variation of phonon energy levels in 
a hybrid system can also affect transitions of photons. So we 
will also study the effect of the ultrastrong coupling between 
the mechanical resonator and the TLS on the photon blockade 
in the hybrid devices. 


III. MASTER EQUATION AND WEAK PUMPING LIMIT 
A. Master equation 


To study photon blockade, we assume that the cavity field 
of the hybrid system is driven by a classical field with fre¬ 
quency U!d, the coupling strength between the driving field 
and the cavity field is |r2|. In the rotating reference frame at 
the frequency iUd, the Hamiltonian of the driven hybrid system 
becomes 


Hr 


hAaa^a + hijJbb^b + —uj^az — hxa'a (b^ + b) 
+hg {b^ +b^ ax + ifi (Ha'^ — H*a) , (8) 


where Aa = uja — uJd describes the detuning between the 
cavity field and the driving filed. 

After introducing the environmental noise, the master equa¬ 
tion of the density operator p for the driven hybrid system can 
be given as 


P — -^[PjHr] + La{p) -b Lb{p) + L„{p). 


(9) 


Here |P) and \M) correspond to quantum states of the TLS, 
and the expression of \P) is given by 


IP) = ^[(l + yi)i/^|t)-(l-v^)^^^U) 


1/2 I 


( 6 ) 


The Lindblad dissipators for the photons and phonons are 
given by 

Eoip) = Xono {opo^ + po — o^op — po^o) 

+ ^ (2opo^ — op — po'^ o) , (10) 


with e = 1 — i^xnjg'f. The expression of \M) can be ob¬ 
tained by replacing | f) and | ],) with | ],) and | f) in Eq. (|6l). 
The state jiy; /3; m) in Eq. © is given as 


where o = a or b corresponds to the variables of the photon 
or phonon, respectively. The Lindblad dissipator for the two- 
level system is 


\ml3;rn) = D{^)S{r])\m). (7) 

The squeezing operator in Eq. 0 is defined as S{r]) = 
exp [^{rjb^'^ — r]*b'^)], while the expression of the displace¬ 
ment operator is D{f]) = exp (/3&i — (3*b). The parameters 
/3 and rj are defined as P{E) = 2ixnEl (hg'^e^, rj = r, and 
exp (2r) = y/e. We can find that the dressed states in the 
hybrid system are formed by the superposition states |P) and 
\M) of the TLS, not the eigenstates | t) or I i) as in com¬ 
mon dressed states. This leads to more complicated phenom¬ 
ena when a probe field passes through such a system. If the 
mechanical resonator and the TLS are in the large detuning 
regime, the energy level spacing between two dressed states 
becomes larger 11^ . 


La{p) = IqTiq {cr_pa+ + a+po^ - a+a^p- pa+a_) 

+ ^{‘2.a-pa+-a+a-p-pa+a-). (11) 

This type of master equations was also studied by Kos- 
sakowski et al. Here 7 a, 76 , and Xq are the decay 

rates of the photon, the phonon, and the TLS, respectively, 
while Ua, nb, and Uq correspond to thermal fluctuation quan¬ 
tum numbers, with rii = l/[exp{huji/{kBT))] (i = a,b,q) 
where ks is, the Boltzmann constant and T is the tempera¬ 
ture. Usually, the thermal photon number Ua can be neglected 
in the low-temperature limit because of the high frequency of 
the cavity held. 

The master equation in Eq. 0 can also be numerically 
solved in the complete basis |n, m, z) (for n, m = 0,1, 2,..., 
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and z =t,Ji) in the case of weak driving field and low tem¬ 
peratures llo 1^ . Because higher excited states can be 
neglected in this case, the photon and phonon numbers can be 
truncated to small values. By numerically solving the mas¬ 
ter equation, we can obtain p which in turn lets us calculate 
various physical properties of the hybrid system. 

B. Analytical solutions in the weak-driving limit 

If the driving field coupling is very weak in comparison 
to the Kerr nonlinearity, and also the temperature is very low, 
then, due to photon blockade, only lower energy levels of the 
cavity field and mechanical resonator are occupied. If the pho¬ 
ton number n and phonon number m are truncated to n = 2 
and m = 1 , respectively, then the quantum state of the hybrid 
system can be written by ll9^100ll 



= Cooi|0,0,;) + Coot|0,0,t) + C'io;|l,04) 

-fCiotll, 0,t) + Coi^lO, 114-) + C'oitlOj Ijt) 
+1^20^12,0,4-) + C'2ot|2,0, t) + Cii4,|l, 1,4-) 

+Ciit|l,l,t) + C2u|2,l,i). (12) 

The coefficients Cnmk (with photon numbers n = 0, 1, 2, 
phonon numbers m = 0,1, and the eigenvalues k =4.,t of 
the dressed TLS states) describe the amplitudes of the corre¬ 
sponding quantum states, andp„j„fc = IC'nmfeP are the corre¬ 
sponding occupation probabilities. 

We use the second-order degree of coherence to describe 
the statistical properties of the cavity field. The equal-time 
second-order degree of coherence is defined by 

In the weak-driving limit, using Eq. (fT2l i and Eq. ( fOl l, the 
second-order degree of coherence can be approximately given 
as 




2 (1(^204.1^ + |C'20tP + 1^^214,1^) 

(|CioiP + |CiotP + |Cn;P + |CntP)2- 


(14) 


The result of Eq. (fT4l i can be used to approximately describe 
the photon statistical properties in the limit of weak driving 
and low temperatures. This will be compared with numerical 
results, calculated using the master equation, in the following 
sections. 

To obtain the coefficients C 104 ,, Cio-t-, Cnj^, Cuf, C 204 ,, 
C' 2 otj C 211 inEq. (fT2l i. we solve the Schrodinger equation 
for the quantum state | ■(/;) of the hybrid system 


.dIV') 

l —;— 

dt 


= HL\ 


(15) 


Here, the effective non-Hermitian Hamiltonian 


= HA'^a^a + Huj'^b^b + — hxa^a + b) 

+hg (fel + b) (Jx + ih (Hal — f2*a) , (16) 


FIG. 3; (Color online) Equal-time second-order degree of coher¬ 
ence as a function of Aa/uJb = (uja — uJd)/oJb in the cav¬ 

ity steady-state limit. The solid curve in each panel is plotted with 
the Rabi model (without RWA), while the dashed curve is plotted 
with the Jaynes-Cummings(J-C) model (with RWA). The curves in 
the three panels correspond to different coupling strengths between 
the mechanical mode and the TLS: (a) gj ( 2 tv ) = 1 MHz; (b) 
gj{2'K) = 5 MHz; and (c) gl{2n) = 15 MHz. The other param¬ 
eters for the three solid curves are: tj(,/(27r) = u}q/{2'K) = 10 MHz, 
7a/(27r) = 0.02 MHz, 76/(27r) = 0.001 MHz, 79/(27r) = 0.002 
MHz, x/(27r) = 0.2 MHz, g/{2'K) = 4 MHz, T = 1 mK, and 
|H|/(27r) = 0.02 MHz. 


includes dissipations with = Aa — Wj = uib — 

and = uig — i'yql2. Here we assume that the ther¬ 
mal fluctuation of the photons, phonons and the TLS can be 
neglected in the extreme low-temperature limit. 

Because we are interested in the statistical properties of the 
cavity field in the steady state, thus we can set d\ip)/dt = 
0. By substituting Eqs. (fT2l i and (fTbl l into Eq. ( fTSl l, we can 
obtain linear equations, as shown in Eqs. dAlTi- dAlOl i of the 
Appendix n By solving these linear equations, we can obtain 
the coefficients in Eq. (fTZt . that is. 


Ciif=inr]2, (17) 

Clot, = C'lO't' = ( 18 ) 

(^204, = (7207 = 12^%: (19) 

C 2 U = (20) 


The expressions of 77^(7 = 1,2,- • -,7) can be found in 
Eqs. (IA.12b . In the weak-driving and low-temperature limit, 
we And that Ciq^, Ciq^, Cnj,, and C\i^ are proportional to fl, 
while C 201 , C 2 o-\-, and C 2 i\. are proportional to The value 
of Coot, will be close to 1 and the amplitudes of the excited 
state tend to 0 if the value of H > 0 . 
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IV. PHOTON BLOCKADE 


In an optomechanical system, the mechanical resonator 
leads to the nonlinearity and energy levels shift of the cavity 
field. For the single-photon state, such shift is Aq = juJb, 
while it is 4An for the two-photo n state 11^ . Both the pho¬ 
ton blockade 11^ and tunneling lIlOlL Il02|] can occur in the 
strong single-photon optomechanical coupling regime. Be¬ 
sides mechanical mode, the TLS can also lead to the variation 
of photon energy levels (see Eq. (|2l)). The energy level struc¬ 
ture of the hybrid system becomes very complex since {ax) is 
a complicated function of x lii. The dressed states formed 
by the TLS and the mechanical mode lead to the splitting of 
phonon energy levels (in standard optomechanical systems). 

We now study how a TLS affects the photon blockade of a 
hybrid system via the second-order degree of coherence 




Tr(pa'f^a^) 
[Tr(pata)]2 ’ 


( 21 ) 


which is calculated here using the master equation in Eq. (|9]l 
and will be compared to the result calculated using Eq. (fT4l) . 
The value of (0) < 1 (0) > 1) corresponds to sub- 

Poisson (or super-Poisson) statistics of the cavity field, which 
is a nonclassical (classical) effect. This effect of the sub- 
Poisson photon statistics is often referred to as photon anti¬ 
bunching. The dips (resonant peaks) of (0) can be used to 
characterize the photon blockade (tunneling) processes. The 
photon blockade describes the single-photon transition, while 
the photon tunneling corresponds to a multi-photon resonant 
transition. 

We plot as a function of Aa/oJb in Eig.[3by using 

the master equation in Eq. (|9]l, the curves in different pan¬ 
els correspond to different values of the coupling strength g 
between the mechanical mode and the TLS. To further study 
the effect of the counter-rotating term on the photon blockade, 
using the master equation in Eq. (|9]), we compare the numeri¬ 
cal results of g^^'> ( 0 ) with (solid curves) and without (dashed 
curves) the RWA. The minimum value of g^^^O) is smaller 
than 1 at the dip near Aa/uib = 0 in the blue solid curve 
of Eig. [3fa), so the photon blockade can be observed. If the 
value of g is much smaller than the transition frequency ujq, 
the TLS has a small effect on the photon blockade in the blue 
solid curve of Eig. Oa) [compare with the black dashed curve 
in Eig. |3c)]. If the value of g becomes larger, the TLS leads 
to two new dips (photon blockade) and several peaks (photon 
tunneling) in the green solid curve of Eig. [Ub). This results 
from the counter-terms -f a-b^) which can be under¬ 
stood by comparing the green solid and brown dashed curves 
of Eig. Ob). When the value of coupling strength g is larger 
than the transition frequency of the TLS, that is p > ujq, more 
dips and peaks appear in the red solid curve of Eig.Oc). And 
the minimum value of g^^'> ( 0 ) near Aa/uib = 0 is larger than 
1 , so the photon blockade in this regime vanishes in ultra¬ 
strong coupling regime. Actually, a similar phenomenon of 
photon blockade in the ultrastrong coupling regime has been 
studied in circuit QED dHt]. 

Eigure O describes the effect of the TLS transition fre¬ 
quency on the photon blockade in the hybrid system. The blue 



-0.1 -0.05 0 0.05 0.1 

A /co, 
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FIG. 4: (Color online) Equal-time second-order degree of coherence 
(0) as a function of Aa/tOb = {oJa — bJd)/ijJb in the steady-state 
limit. The curves in the three panels correspond to different transition 
frequencies of the TLS: (a) ojq/ (27r) = 10 MHz; (b) tOq/ (27r) = 13 
MHz; and (c) LUq/ {2 tt) = 20 MHz. The other parameters are as in 
Fig.m except gl{2Tv) = 4 MHz. 


solid curve of Fig. Ufa) describes g^^^{0) when the mechani¬ 
cal resonator interacts resonantly with a TLS, while the green 
solid [in Figs. HJb)] and red solid curves [in Fig. HJc)] corre¬ 
spond to the detuning cases. When the mechanical mode and 
the TLS are in the detuning regime, the positions of the left 
and right dips (relative to the point Aa/uib = 0 ) are changed 
in the green solid curve of Fig. lUb). The minimum value of 
g^^^ ( 0 ) of the left dip becomes larger than 1 , so the photon 
blockade disappears near this point. But the photon blockade 
near the right dip is enhanced. If the detuning \uiq — ujb\ is 
larger than the coupling strength g, all the dips and peaks in¬ 
duced by the TLS disappear in the red solid curve of Fig.Hfc), 
in this case the photon blockade is similar to that of stan¬ 
dard optomechanical systems [see the black solid curve in 
Fig. Sc)]. 

The effect of the decay rate 7 , on the photon blockade of 
optomechanical systems is discussed in Fig. If the value 
of 7 q becomes larger, the dip near Aa/uib = 0 is almost 
invariant, but the left and right dips (relative to the point 
Aa/wfc = 0) change greatly. The minimum value of g^^\0) 
near the right dip is even larger than 1 , so the photon blockade 
disappears in this regime. The photon blockade near the left 
dip will also vanish if the value of 7 , continues to increase. If 
the decay rate becomes very large, all the new dips and peaks 
induced by the TLS vanish, and the photon blockade in the 
red solid curve of Fig.|5lc) is then almost then same to that of 
standard optomechanical system [see the black dashed curve 
in Fig. He)]. 

In Fig. | 6 l we compare the results obtained numerically by 
the master equation in Eq. (|9]) with those obtained analyti¬ 
cally in Eq. ([T4l) . We plot as a function of the detun¬ 

ing Aa/ojb for r = 0 K. The blue solid curves in Figs.| 6 l[a) 
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FIG. 5; (Color online) Equal-time second-order degree of coherence 
(0) as a function of Aa/ujb = (uJa — ujd)li-Ub in the steady-state 
limit. The curves in three panels correspond to different decay rates 
of TLS: (a) 79/(27r) = 0.002 MHz; (b) 'yg/{2n) = 0.005 MHz; 
and (c) 7q/(27r) =0.1 MHz. Here gj{2TT) = 4 MHz for the solid 
curves, while gl{2n) = 0 Hz for black dashed curve in panel (c), 
and the other parameters are same as in Fig.[3 


and |6|b) are plotted using the master equation, while the red 
dashed curves are plotted with the analytical result in Eq. (fT4l) . 
The results of two methods are almost the same for Fig.|6la). 
If the coupling strength g becomes larger, the deviations be¬ 
tween the blue solid and red dashed curves in Fig. |6|b) be¬ 
comes larger. This difference originates from the approxi¬ 
mation when Eq. (fT4l) was derived, because, contrary to our 
precise numerical calculations, we have assumed in our ana¬ 
lytical approach that (i) the system evolution is pure and (ii) 
some transition processes, such as |0,2,4,), |0,2,t), |1,2,4,), 
etc., were neglected. 

Therefore, we conclude that the coupling strength (to the 
mechanical mode), transition frequency, and the decay rate of 
the TLS can be used to tune the photon blockade and tunneling 
of optomechanical systems. 


V. CONCLUSIONS 

We have studied single-photon blockade and tunnel- 
ing(corresponding to multi-photon blockade) of a hybrid sys¬ 
tem consisting of an optomechanical cavity and a TLS. We 
find that the photon blockade of the optomechanical device is 
significantly affected by a TLS when it is coupled to the me¬ 
chanical resonator. Compared with the results of only the op¬ 
tomechanical cavity, the TLS shifts and splits the peaks and 
dips of the second-order degree of coherence of the cavity 
field in the optomechanical subset. We also find that the TLS 
gives rise to several new peaks and dips in the second-order 
degree of coherence of the cavity field in the hybrid system. 

If the coupling strength (between the mechanical mode and 




FIG. 6: (Color online) Equal-time second-order degree of coherence 
g^'^\0) as a function of Aa/oJb at zero temperature in the steady- 
state limit. The blue solid curves correspond to our numerical precise 
solutions of the master equation, while the red dashed curves are 
plotted with the analytical result in Eq. 01- Both figures are plotted 
by substituting with a+(T- in Eq. (1^. Here x/{2'ir) = 0.5 MHz, 
|fl|/(27r) = 0.01 MHz, and T = 0 K, while other parameters are 
the same as in Eig.[^ 


the TLS) is comparable or larger than the transition frequency 
of the TLS, new blockade dips and resonant peaks appear for 
the second degree of coherence. Moreover the new blockade 
dips and resonant peaks can be tuned if we change the transi¬ 
tion frequency or the decay rate of the TLS. The photon anti¬ 
bunching of hybrid systems can also be tuned if we change 
the parameters of the TLS. That is, our study may provide a 
new method to control and tune the nonlinearity and nonclas- 
sical effect of the cavity field of the optomechanical system by 
coupling to a tunable TLS. Our calculation may also provide 
an approach to detect a low-frequency TLS (which might be a 
defect in a mechanical resonator) using optomechanics. 
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The equation 0 = —w'( 7004,72 — tO*( 7 io 4 , has no physical 
meaning if |0| —0, so it can be neglected. The system 
has a probability to remain in the ground state | 0 , 0 , 4 ,) in the 
weak-pumping limit, so we can set ( 7 oo 4 , = 1 (then the expan¬ 
sion coefficients are unnormalized). The terms (zv/ 2 O*( 7204 .) 
in Eq. (IA.4b . (iO*( 72 ot) in Eq. (IA.51 l. and (—z\/ 20 *( 72 i^) in 
Eq. ( IA. 6 I 1 are of higher-order in 0, so they can be neglected. 

Through some calculations, we can obtain the solutions of 
the expansion coefficients in Eq. (fT2T i. The corresponding co¬ 
efficients rii{i = 1,2,... ,7) are given by 

m = (xAT’Aii- f pAT’Ag - A^pAg)/£>i, 

V 2 = — (^A^Ag-b (7A2 As + /-Di, 

m = (x7i -572 - 1 )/A^, 

m = (x72 - A4Pi)/A3, 

V5 = (A13A14-b 2 xpAi 5 )/ZJ2, 


Appendix: The expansion coefficients of Eq. 

In the weak pumping limit, only low excited states of pho¬ 
tons and phonons are occupied. If the temperature is ex¬ 
tremely low, the quantum state of the hybrid system can be 
written as a sum of the finite orthogonal basis states given in 
Eq. (IT 2 I 1 . Considering the effect of environment noises, the ef¬ 
fective non-Hermitian Hamiltonian of the hybrid system can 
be obtained in Eq. (fTfil l. In the steady state case, we can set 
d\ijj)/dt = 0. Substituting Eqs. (fT2T i and (fTfil l into Eq. (fTSl l. 


rje — — (A12A15 -b 2 xpAi 4 ) / D2, 

dr = 2% (A13A14-b 2xpAi5)/(AgUg) + x/^pi/Ag 

(A12A15-b 2x3X14)/(A5U2), (A. 12 ) 


with Di = A 7 A 2 (AgAio — AgAii) and D 2 = A 12 A 13 — 
4x^3^. The parameters Xi{i = 1,2,..., 15) in the above 
equations are 


Eq. (fT^ as follow: 


Ai 

= Ar - 232 /^;, 

0 = (wq/ 2 ) (7oot + ffC' 014 , — iH*(7io-t-, 

(A.l) 

A 2 

= A 2 + - 2|H| 

0 = A7(7oi4. -b 3(7oot ~ *^*C'ii4,, 

(A.2) 

A 3 

= A 2 - 4|0| V/ (Aiu;;^) , 

0 = AT'Coi't-+ 3(7 oo4, — *ll*(7ii't-, 

(A.3) 

A 4 

= 3 -b 23 Hp/ (Aiw') , 

0 = A 2 (7io 4, — xCiij, + 3C'iit + *^^* 004 . 


As 

= A 3 - - |H|VAi, 

-iV2n*C20i, 

(A.4) 

Ae 

= 3 -b 2 U ^ 3 /Aia;', 

0 = A^(7io-|'— xC'iit + '?C'ii 4 . + iO(7Qo^ 


A 7 

= A 3 + - | 0 |VA+, 

-iV2n*C20t,, 

(A.5) 

As 

= A 5 — X^/A 2 — A 4 A 6 /A 3 , 

0 = A;j'(7ii4, — xCio4, + ^Ciot + *A1(7oi4. 


Ag 

= 9X/^2 TxAe/Ag, 

-iV2n*C2n,, 

(A. 6 ) 

Aio 

= 9X/^2 +XA 4 /A 3 , 

0 = A;7C'ii-|' ~ xC'iot + <?C'io4. + iH(7oit5 

(A.7) 

All 

= A 7 — x^/Ag — 3^/A2 , 

0 = A 4 (7204. — 2x(7214, + t'\/2HCio4,, 

(A. 8 ) 

Ai 2 

= A4A5-4x^ 

0 = A 4 (72ot + 9^211 + i-\/2H(7io-t-, 

(A.9) 

Ai3 

= A+A 5 - 3 ", 

0 = A5(7214, — 2x( 7204, + 3(72ot + iV2rtCu 

i^A.lO) 

Ai4 

= \/2(2x7i + A 533 ), 

with the definitions of Ai{i = 1, 2,..., 9) being 


Ais 

= V2grii - v/2A534. 


(A.13) 


T 

T W^/2, 

T ^^g/2, 

AJ = 2A^ W5/2, 

^5 = 2A7 + wTi - w'/2. 


With the expressions of 31(1 = 1, 2,..., 7) in Eqs. (IA.121 i. 
we can obtain the analytical result of the second-order de¬ 
gree of coherence in Eq. (fl4l i in the weak-pumping and low- 
temperature limit. 
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